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g_) , Abstract 

We use semiclassical methods to calculate the probability of induc- 
ing a change of topology via a high-energy collision in the S'C/(2)-Higgs 
^ ■ theory. This probability is determined by a complex solution to a clas- 

sical boundary value problem on a contour in the complex time plane. 
In the case of small particle number it is the probability of instanton- 
like processes in particle collisions. We obtain numerically configu- 
rations with the correct topological features and expected properties 
in the complex time plane. Our work demonstrates the feasibility of 
the numerical approach to the calculation of instanton-like processes 
in gauge theories. We present our preliminary results for the suppres- 
sion factor of topology changing processes, which cover a wide range 
of incoming particle numbers and energies below the sphaleron energy. 
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1 Introduction and Summary 



Non-perturbative effects occur in many processes studied by quantum field 
theory. Well known examples are the decay of the false vacuum in scalar 
models and instanton-like transitions in gauge theories. The latter 
companied by non-conservation of fermion quantum numbers 0, such as 
baryon number, and thus are of importance for particle physics and cosmol- 
ogy. 

In weakly coupled theories, these processes are described at low energies 
by classical Euclidean solutions interpolating between initial and final states 
separated by an energy barrier. In the examples mentioned above such solu- 
tions are known as bounce and instanton 0. The Euclidean action of the 
solution determines the exponential part of the transition rate. It is inversely 
proportional to a small coupling constant present in the model and thus, in 
general, the rate is highly suppressed. 

At energies comparable to the height of the barrier, the probabilities of the 
transitions between topologically distinct vacua may become unsuppressed. 
This takes place, for example, at finite temperature 0, ^ finite fermion 
density J^, ^, or in the presence of heavy fermions in the initial state 
PH] , [TT|, |r^. In high energy particle collisions the situation is not quite the 
same. 



As was first noted in |T3|, |14[, at relatively low energies the corrections to 
the tunneling rate can be calculated by perturbative expansion in the back- 
ground of the instanton (bounce). Further studies showed that the actual 
expansion parameter is e = E/Egph and the total cross section 

of induced tunneling has an exponential form 

atot{E) ~ exp 1-^Fhg{E/E,p,,) 

where g is the small coupling constant and the function Fhg{^) is a series in 



powers oi e = E / E^^^^ (for a review see |[T9||). 

While the perturbation theory in e is limited to small e, the general form 
of the total cross section implies that there might exist a semiclassical-type 
procedure which would allow, at least in principle, to calculate Fhg{£) at 
£ > 1. Since the initial state of two highly energetic particles is not semi- 
classical, the standard semiclassical procedure does not apply and a proper 
generalization is needed, which was proposed in refs. [EO, pTl [2^. The cor- 



responding formalism reduces the calculation of the exponential suppression 
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factor to a certain classical boundary value problem, whose analytical solu- 
tion is not usually possible. 

The semiclassical approach proposed in refs. |21|, is based on the 
conjecture that, with exponential accuracy, the two-particle initial state can 
be substituted by a multiparticle one provided that the number of particles 
is not parametrically large. The few-particle initial state, in turn, can be 
considered as a limiting case of truly multiparticle one with the number of 
particles = v/g^ when the parameter v is sent to zero. For the multi- 
particle initial state the transition rate is explicitly semiclassical and has the 
form 

AT) ~ exp I — ^F{e,u) 

According to the above conjecture, the function Fnci^), corresponding to 
the two-particle incoming state, is reproduced in the limit z/ — 0, 

hmF{e,u)=FHGie). 

Therefore, although indirectly, the function Fhg{.^) is also calculable semi- 
classically. Although not proven rigorously, this conjecture was checked ex- 
plicitly in field theory models in several orders of perturbation theory in e PB 
and also in a quantum-mechanical model for all e 



Until now the only analyses of induced tunneling have been performed 
in a quantum mechanical model and in a scalar field model of false vac- 



uum decay A particularly interesting case to study is, however, the 

Electroweak Theory where different topological sectors are separated by a 
potential barrier of the height i?sph ~ lOTeV [^, Whether the exponen- 
tial suppression disappears in this theory at sufficiently high energy is still 
an open question. In this paper we consider an SU{2) gauge theory with 
Higgs doublet, which corresponds to the Electroweak sector of the Standard 
Model with Ow = 0. 

Classically allowed over-barrier sphaleron transitions in this model were 



studied in ref. [^. All solutions found in ref. |28] are configurations with 
large number of particles in the initial state and thus they do not correspond 
to few-particle collisions. 



In the present work we adapt the prescription of |2y, 0, ^ to theories 
with gauge degrees of freedom. The prescription requires the solution of field 
equations on a contour in complex time plane (see figure |lD with boundary 
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conditions of a special form. Since the solutions cannot be found analytically, 
one has to invoke numerical techniques. Here we present solutions to this 
problem in a limited region of parameters. Our solutions possess all expected 
features, including correct topological structure and singularity structure in 
complex time plane. Thus, our results appear to validate the use of numerical 
methods for the study of instanton-like transitions in gauge theories at high 
energies. 

In this work we explore the region of parameters with E < Esph and 
0.4 <!/<!. We calculate the exponent of the suppression factors for a 
wide range of values of E and u within the above region. Not surprisingly, 
since we work at energies below the sphaleron energy, we find that topology 
changing processes remain exponentially suppressed throughout the region 
we studied. More relevant may be the fact that the trend of our results 
appears to indicate that topology changing processes with low initial particle 
number may remain suppressed well above -Egph- 

The paper is organized as follows. In section ^ we briefly review the 
boundary value problem for the semiclassical calculation of a{E,N). In 
section ^ we present the SU (2) model used in our calculations. In section 
^ we describe our numerical results and section ^ contains a few concluding 
remarks. 



2 Semiclassical approach to induced tunnel- 
ing probability 

The inclusive multiparticle probability of a transition from a state with fixed 
energy E and number of particles about one vacuum to any state about 
another vacuum can be written in the form: 

a{E,N) = J2\{f\SPEpN\t)\\ 

where 5" is the 5'-matrix, Pe,n are projectors onto subspaces of fixed energy 
E and fixed number of particles A^, and the states \i) and |/) are perturbative 
excitations about topologically distinct vacua. The method of semiclassical 



calculation of this probability a{E, N) was formulated in refs. pO| , |21| , |22| , |25 
Here we only review the prescription. 

For small coupling constant g, the semiclassical approximation is appli- 
cable. It boils down to the classical boundary value problem specified on the 
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Figure 1: The contour in complex time plane used in the formulation of the 
boundary value problem (^) . Crossed circles represent singularities of the field at 
r = 0, for larger r the singularities generally move to larger times. 



contour in complex time plane shown in figure |I]. In the initial part of the 
contour (AB) at sufficiently large negative times the fields have to be in the 
linear regime. Let us denote the fields collectively as (p{x,t). The frequency 
components at distant past (part A of the contour), and g,,, are defined 
as follows, 

'^^'^'^^l^^-y (2vr)3/2^W +9ke ) 

The tunneling probability reads (after proper rescaling of the fields) 

aiE,N)r.exp!^-^F{s,u)^ (1) 
-F{e, u)=NQ + ET - 2lm[SABCD{v)] +ReBi , (2) 



9' 



where 



(it is easy to check that the expression (H) for F{e, v) is independent of Tj if 
the system is in linear regime at initial time). Here the field (p interpolates 
between neighborhoods of topologically distinct vacua and satisfies the field 
equation: 
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At initial time the frequency components of the solution should satisfy the 
following equations ( "0 boundary condition" ) 



h = ^-""dk ■ (3b) 

On the final part of the contour (CD), the field is real, so that 

Im0(a;,O) = lm(^(a;,0) = . (3c) 



Equations (paD-(3c) specify the boundary value problem corresponding to 
the induced topological transition. 

The quantities entering equation (g) are defined as follows. Sabcd is 
the action for the solution of equations (|^), energy and number of incoming 
particles are 

E = Jdk uj^gl , N = Jdk Ugl . (4) 

These equations indirectly fix the values of the auxiliary variables T and 
for given energy and number of particles. Alternatively, one can fix T and 0, 
solve the boundary value problem (^ and obtain the corresponding values of 
E and using (^. This is especially convenient in numerical calculations. 

The interpretation of solutions to the boundary value problem (^ is as 
follows. On the part CD of the contour the saddle-point field is real; it 
describes the evolution of the system after tunneling. On the contrary, it 
follows from boundary conditions ([3^1) that in the initial asymptotic region, 
the saddle-point field is complex provided that 7^ 0. Thus, the initial state 
which maximizes the probability (^ is not described by a real classical field, 
i.e. the classical field must be analytically continued to complex values and 
this stage of the evolution is essentially quantum even at A^ ~ 1/(7^. 

The picture described implies that there exist singular points of the so- 
lution in the complex time plane, as shown in figure 0. To understand this, 
one notices that on the negative part of the real time axis, the solution (at 
least for energies below the sphaleron energy) "bounces back" to the same 
vacuum, as in the CD part. On the other hand, the solution in the AB part 
of the contour and its analytic continuation to the real axis is close to a dif- 
ferent vacuum. This may happen only if a branch cut exists between the real 
axis and AB part of the contour. Similar arguments require a singularity to 
exist to the right of the BC part of the contour too. 
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The case = is exceptional. In this case, the boundary condition ( ^BD 
reduces to the reahty condition imposed at Imt = T/2. The solution to 
the resulting boundary value problem is a periodic instanton of ref . |^ . The 
periodic instanton is a real periodic solution to the Euclidean field equations 
with period T and two turning points at t = and t = iT/2. Being ana- 
lytically continued to the Minkowskian domain through the turning points, 
periodic instanton is real at the lines Imt = and Imt = T/2 and therefore 
satisfies the boundary value problem (0) with B = 0. Periodic instanton 
solutions have been studied with a computational approach similar to the 
one used in this paper in ref. [pO| . 



3 The model 

In this paper we study the four-dimensional model which captures all the 
important features of the Standard Model — an SU{2) gauge theory with the 
Higgs doublet. This model corresponds to the bosonic sector of the Standard 
Model with 9y/ = 0. Also, according to ref. we ignore the back reaction 



of fermions on the gauge and Higgs fields dynamics. The action of the model 
is 

S = y"d^x|-iTrF^,F^'^ + (D^<l>)tD^$-A($t$-l)2| , (5) 

where 

D^^ = {d,-iA^)^ 

with A^ = A1^a"'/2. Here we have eliminated some inessential constants by 
an appropriate choice of units. We have also set the gauge coupling constant 
g = 1 hj proper rescaling of the fields and action, but it can be easily 
restored in the final result (see (0)). The Higgs self-coupling A was set equal 
to A = 0.125 in all calculations, which corresponds to ttih = mw 

This theory is still too complicated for numerical study because of large 
number of variables. To make the problem computationally manageable 
we consider only configurations spherically symmetric in space |3^, which 
reduces the system to an effective 2-dimensional theory. It still possesses 
many features of the full 4-dimensional model, such as a similar topological 
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structure. Moreover, for large times the energy disperses along the radial 
direction as it would do in the full 4-dimensional theory. This guarantees 
that the system linearizes with time making it possible to impose boundary 
conditions (pB]) in the asymptotic region. 

The spherical Ansatz is given by expressing the fields in terms of six real 
functions Oq, cti, a, ^ and u: 



Ao{x,t) = -ao{r,t)a- • n 



, . a(r,t) 
ai{r,t)(T • nui H (a^ 



1 + P{r,t) 



(6 



b) 



<l>(a;, t) = [/i(r, t) + iu{r, t)(T • n]^ 



(6c) 



where n is the unit three-vector in the radial direction and is an arbitrary 
constant two-component complex unit vector. The action (|^) expressed in 
terms of the new fields becomes 



S = 47r dt 



dr 



1 



XX 



if-lixx + i 



where the indices /z, v run from to 1 and 



X = a + ip 
(p = fi + iu 
Df,x = (df, - ia^)x 



X = a-i(3 
(j) = fi — iu 
Df,x = (df, + ia^)x 

— — % — 



(7a) 
(7b) 
(7c) 
(7d) 

(7e) 



Note that the overbar on (p, x and denotes changing z — z in the 
definitions (|^) above, which is the same as complex conjugation only if the 
six fields a^, a, P, ^ and z/ are real. In the boundary value problem (|^) these 
fields become complex and overbar no longer corresponds to normal complex 
conjugation. 
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Vacuum structure. The spherical Ansatz (|^) has a residual U{1) gauge 
invariance 



a^t ^ ct/. + d^Q. 
X - 

with gauge function VL{r,t). The complex "scalar" fields x ^-nd have U{1) 
charges 1 and 1/2 respectively, is the f/(l) gauge field, f^i, is the field 
strength, and in (|^) is the covariant derivative. The residual U{1) gauge 
invariance must be fixed when solving the equations numerically. In our work 
we chose the temporal gauge oq = and impose Gauss' law (the equation 
corresponding to variation over oq) at the initial moment of time. 
The trivial space-independent vacuum of the model has the form 

Xvac ^ ; 0vac il ; '^1 vac . 

Other vacua can be obtained from the trivial one by the gauge transforma- 
tions: 

Ct/ivac dfj^a 
Xvac ^6 
Vvac -"-"^ 

Q should be zero at origin. Vacua with different winding numbers correspond 
to — 2n7r as r — >■ oo. For such values of Q, the fields of the original four- 
dimensional model are constant at spatial infinity, which is the standard 
choice. It allows for a standard description of the topological properties 
of vacua — since the sphere S"^ at spatial infinity is mapped to one point 
in field space, one can compactify the space to 5*^ and consider mappings 

— >• SU{2) (the latter correspond to pure gauge field configurations). 

One can also make any other choice of fields at spatial infinity (as long 
as the fields are pure gauge and constant in time there). In our case it is 
convenient to set fl — *■ {2n — l)7r at r ^ oo. This is equivalent to the 
requirement that the fields satisfy the following boundary conditions at the 
"boundaries" of space: 

X|i — >0 +00 ^ i 

RedrCplr^O 0|r-*oo — > i 

Im0L,^n . 
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r = oo 




(a) (b) (c) 

Figure 2: Topological transition in the SU{2) Higgs model: behavior of the fields 
(j) and X- Bold arrows show the change of the field as the radial coordinate increases 
from r = to r = oo. The configurations are shown: (a) at initial time, (b) in the 
middle of the process and (c) at final time. 

In the original 4- dimensional theory this means that the sphere S"^ at spatial 
infinity is mapped onto the equatorial sphere of SU{2) parameterizing the 
gauge vacua. 

In this gauge no r-independent vacuum exists, but transition from vacua 
with n = and = 1 is described in a very symmetric way. The behavior 
of the fields x ^^id (p for such transition is shown in fig. In the original 4- 
dimensional model this topology changing process corresponds to a transition 
where the fields wind over the lower hemisphere of SU (2) before the transition 
and over the upper hemisphere after the transition. 

Boundary conditions. The formulation of the G boundary conditions 
([3bD requires additional effort in this model, as compared to the case of single 
scalar field. The reason is that one of the fields is not really a physical field 
in = gauge — if fixed at some time, it can be expressed in terms of other 
fields at any other moment of time using Gauss' law, which is a first order 
equation (in time). The corresponding expression is quite complicated, so it 
is comfortable instead to solve the second order equations only and impose 
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Gauss' law at one moment of time. Then the Gauss' law is automatically 
satisfied at all other times. 

So, we have to impose boundary conditions ( pE| ) only on four of the 
five fields ai, a, (3, fi, v (or, to be more precise, on four combinations of these 
fields). At one moment of time we also have to impose Gauss' law constraint 
and fix the gauge freedom to make the boundary value problem properly 
defined. The two latter constraints are equations on complex functions of 
spatial coordinate r (recall that all fields in our problem are complex and 
gauge function may be complex too). Half of these equations are not 
needed — they are duplicated by the reality conditions at the CD part of the 
contour. Reality implies that imaginary part of Gauss' law is zero and forbids 
gauge transformations with imaginary fi. So we are left only with the real 
part of Gauss' law and have to fix invariance under real gauge transformations 
at the initial time moment. The latter can be done by setting a certain 
combination of fields corresponding to unphysical perturbation of the initial 
vacuum to zero. This gives us the needed number of boundary conditions to 
determine unique (in general) solution. 

Zero mode. One more complication is that, in continuous formulation, the 
boundary value problem (j^) does have an invariance under translations along 
real time (both field equations and the boundary conditions are invariant 
under such translation). To properly define the boundary value problem one 
has to fix the position of the solution in time. In the lattice version this 
invariance is violated by the discretization and finite volume effects. 

A possible modification of the equations is the following. One of the 
equations (^) (in lattice case k takes discrete values) is changed to 

I4l = e-®|^7fc|. 

If the field is in linear regime at initial time relative phase between and 
is zero because total energy is real. So in linear regime we get the solution 
to the original boundary value problem. 

Instead of the equation for the relative phase any (real) equation which 
is not invariant under time translations can be used. It fixes the position of 
the solution. We fix the center of spatial distribution of field x at part A of 
the contour, which corresponds to fixing of the position of incoming wave. 
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4 Numerical results 



There are several peculiar features of the boundary value problem (^) that 
make solving it numerically a computational challenge. First, equation (|3al) 
is non-linear. Second, at O 7^ the fields are necessarily complex and the 
solution is not a maximum of (|^) but only a saddle point. Third, the time 
contour contains both Minkowskian (AB and CD) and Euclidean (BC) parts, 
so the problem is both of hyperbolic and elliptic type. Finally, the initial 
boundary conditions (|3BD should be imposed when the system is already 
close to linearity which requires large spatial volume of the configuration. 

These factors constrain lattice parameters and applicable numerical tech- 
niques. In the lattice version the boundary value problem (^) becomes a set of 
non-hnear algebraic equations for the field values ipij = {ai, a, (3, fi, i'}(ti, Vj) 
at the lattice sites with coordinates (ti,rj), where tq, . . . ,rN corresponds to 
spatial radial direction, . . . yt^t+i are complex time coordinates on con- 
tour ABCD, index J is a combination of the spatial index j and field type 
{ai, a, /3, /X, I/} and runs from 1 to 5A^ — 4. The lattice size is characterized 
by the lengths of AB, BC and CD parts of the time contour (T/vfinitiai, Te 
and TMfinai respectively), and spatial size L. While on the Euclidean part 
the solution is compact in space (and has characteristic size of order of 1), it 
is generally evolving along the light cone in the Minkowskian regions. This 
requires that L > Tjv/initai- In calculations we chose L = 8 and T^/initai = 6 
(in units of Higgs boson mass). 

The numerical method to solve the set of equations which constitute the 
lattice version of the boundary value problem (|^) is to be chosen according to 
the specifics of the problem described above. To get rid of the non-linearity 
we employ a multidimensional Newton-Raphson method which approaches 
the desired solution iteratively. At each iteration, the linearized equations 
in the background of the current approximation are solved. The next ap- 
proximation is obtained by adding the solution to the background, and the 
procedure is repeated. The advantage of the algorithm is that it does not 
require positive-definiteness of the matrix of second derivatives. It is, how- 
ever, sensitive to zero modes. In the absence of zero modes, the algorithm 
converges quadratically; the accuracy of 10~^ is typically reached in 3-5 it- 
erations. The convergence slows down in the presence of very soft modes, as 
typically happens near bifurcation points. 

At each Newton-Raphson iteration one solves the set of Nt x 5A^ linear 
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equations of the general form 



L ■ u = d 



where u is the vector formed of A^^ ■ 5A^ unknowns, L is the matrix of dimen- 
sion Nt ■ 5N X Nt ■ 5N (first variation of the full non-linear equations) and d 
is a constant vector (full equations evaluated at the current background; at 
the desired solution d = 0). The inversion of this matrix is the most time 
consuming part of the calculation; its efficiency determines how large Nt and 
can be used. The matrix L is neither positive-definite nor even symmet- 
ric, but has a special sparse structure as it originates from the second order 
differential equation. The linear equations relate only adjacent time slices 
i — i, i + so one can eliminate equations for alternate time slices ^4, 30 



Moreover, one can do that in parallel, making use of the power of multipro- 
cessor computers. The algorithm requires ~ Nt{5N)^ multiplications. Note, 
that it is highly asymmetric in Nt and N — one can use large A''^ but is very 
constrained with the choice of A^. The calculations were performed with 
A^ = 64 and A^"* = 350. 

The Newton-Raphson method requires a good initial approximation for 
the solution. This favors the following general strategy. We first find the 
periodic instanton solution (which corresponds to = 0) After the 

periodic instanton is found, we change parameters T and by small steps, 
using the solution from the previous run as a starting configuration for the 
next one. At each step we calculate the energy E, number of particles A^ and 
the exponential suppression factor F(e, i/). 

A typical configuration is shown in figure |^. One can clearly see that the 
phase of the field x behaves as shown in figure going from — vr at r = to tt 
at r = L along different sides of the circle in the initial and final states. In the 
middle of the Euclidean region, a zero of the field x is present (center of the 
"instanton"). Also the incoming/out going wave is observed at initial/final 
Minkowskian part. 

It is fairly straightforward to check that the singularities shown in fig- 
ure (|I|) are indeed present, by continuing the solution from the part BC of 
the contour to the whole complex time plane. The left singularity is approx- 
imately at the same distance from the Euclidean part of the contour for all 
solutions. The right singularity moves to larger positive times as the energy 
increases. 

Summary of the results is given in figure ^ Lines of constant suppression 



13 



Figure 3: Configuration for r/2 = 1.98 and 9 = 2.74. e = 0.711. Heiglit and 
color represent |x| and the phase of X; respectively. For visualization purposes the 
Euclidean part of the contour is inclined. 



factor F are shown on e-z/ plane. Units are normalized to the sphaleron: 

Behavior of the constant suppression lines has the following features. Near 
the periodic instanton (G = 0), the dependence of F on z/ is weak. This can 
be seen analytically from the boundary problem (|^) itself. Making use of the 
fact that F is stationary with respect to T and B one finds 



dN 
dE 



F=const 



T 

e 



which is infinite as G — 0. 

When one moves away from the periodic instantons, lines of constant F 
flatten out; in other words, increase of energy in this region leads to smaller 
decrease of the suppression exponent than in the vicinity of the periodic 
instanton. 

The quantity of primary interest is the two-particle cross section 
atot(F) ~ exp I — —F {E/E,^^, 0) 
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Figure 4: Lines of F{e, v) = const. F is normalized to be unity at = e = 
(instanton). Diagonal line directed from sphaleron {v = e = 1) towards the origin 
is the line of periodic instantons. Open symbols correspond to configurations with 
wrong topology. Numbers near the curves show the suppression factor F. 



We can plot z/) as a function of z/ for different values of energy (fig. ^. 
Extrapolation of the data for e = 0.4, . . . , 1.0 crosses the line u = at 
F ~ 0.8, . . . , 1. So we conclude that for energies below the sphaleron, the 
exponent of the suppression factor of the tunneling processes is only about 
20% less than for zero energy (instanton case). 

Finally, the open symbols in fig. ^ correspond to evolutions where the 
field configuration, after having undergone a topological transition (marked 
by the zero of the x fi^ld on BC part of the time contour — see fig. |^), goes 
trough a reverse transition on the CD part of the contour and returns to 
the original topological sector. We will discuss these solutions in the next 
section. 
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Figure 5: Dependence of the suppression factor F on for different energies. 
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5 Conclusions and Outlook 

Of course, the point of major interest is whether at sufficiently high energy 
particle collisions can induce unsuppressed topological transitions and baryon 
number violation in the Electroweak Theory. In terms of the graph of fig. | 
the question is whether the line corresponding to F{e, z/) = does approach 
the t/ = axis asymptotically for increasing e and, if it does, at what rate. 

It should be noted that in our study we have not yet been able to obtain 
topology changing solutions with F = 0. As we described, the exploration of 
the space of solutions is done through a deformation procedure, by which the 
parameters T and 6 are gradually changed. This in turn moves the solutions 
in the e, v plane. We observed that as one goes beyond the sphaleron energy, 
the solutions become unstable. The field configurations linger for longer and 
longer times in the neighborhood of the sphaleron and then return to the orig- 
inal topological sector. One thus obtains solutions satisfying the equations 
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of motion and all the boundary conditions, but for the requirement that the 
topology changes. This instability is not unexpected. It was observed by the 
authors of ref. |^ in their study of transitions across a potential barrier in a 
quantum mechanical model. It indicates a likely bifurcation in the space of 
semiclassical solutions as one reaches the barrier energy. In ref. [^] the prob- 
lem was solved by deforming the time contour for the evolution (see fig. |lD 
to make it go below the real axis (the segment CD in fig. |l|) before returning 
to the real time axis. This helped pinpoint a set of solutions with transition 
across the barrier. We are pursuing a similar strategy in our current search. 
The computational problems one faces are nevertheless daunting. The most 
crucial factor is the ability of following the solutions for negative and positive 
time well into the linear regime, where they are settled in the two different 
different topological sectors. This in turn requires the use of a grid of large 
extent in the radial direction and makes the computation quite demanding. 
We are making progress and hope that we will be able to report on topol- 
ogy changing solutions above the sphaleron energy in the near future. For 
the moment, we believe that the detailed information we obtained for the 
lines of constant suppression factor below the sphaleron energy can also be 
of value. In particular, the marked bending of the lines toward increasing 
energy as one lowers the incoming particle number seems to indicate that 
topology changing transitions in particle collisions will occur, if at all, only 
for energy substantially higher than -Egph- 
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